Quantum-classical correspondence of entropy contours in the transition to chaos 



Raphael Zarum and Sarben Sarkar 
Department of Physics, King 's College London, Strand, London WC2R 2LS, UK 

(February 5, 2008) 

Von Neumann entropy production rates of the quantised kicked rotor interacting with an en- 
vironment are calculated. A significant correspondence is found between the entropy contours of 
the classical and quantised systems. This is a quantitative tool for describing quantum-classical 
correspondence in the transition to chaos. 
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I. INTRODUCTION 



The fundamental split between integrable and non- 
integrable systems in classical mechanics has not been 
comprehensively mirrored in quantum mechanics |Q. 
The issue seems to hinge on finding a suitable definition 
for quantum chaos. The sensitive dependence to initial 
conditions that characterises classical chaos is wholly un- 
derstood in terms of trajectories of classical phase space 
points which have no direct quantum analog H . This has 
led to two distinct ways of identifying variables to mea- 
sure quantum chaos. One involves investigating various 
quantum variables that can act as signatures of chaos by 
clearly distinguishing between quantum systems whose 
classical counterparts are integrable and those which are 
non-integrable. The implementation of an expanding ar- 
ray of energy spectra properties have made this approach 
highly successful (e.g. J^-fj]])- A second approach seeks 
an intrinsically quantum definition of quantum chaos by 
investigating the quantum parallels for variables such 
as Lyapunov exponents and various entropy measures 
which define and quantify chaos in classical mechanics 

BUS 

In this paper we adopt the second approach and de- 
velop an original technique involving the analysis of en- 
tropy production measures to reveal a clear correspon- 
dence between the quantum and classical formulations of 
a seminal system which has a rich structure of delicately 
interwoven regular and chaotic dynamics - the standard 
map. Section II A briefly outlines the main characteris- 
tics of the standard map and then Section II B describes 
how classical entropy contours are calculated and used 
to give a comprehensive account of these characteristics. 
In a similar fashion, Section III A briefly outlines the 
quantisation of the standard map and then Section III B 
describes how quantum entropy contours are generated 
through interaction with an environment. The similar- 
ities and differences between the classical and quantum 
entropy contours are explained in Section IV A and Sec- 
tion IV B concludes the paper with a discussion on the 
use of entropy measures to describe quantum chaos. 



II. CLASSICAL STANDARD MAP 



A. Phase space description 

The standard map describes the local behaviour of 
nonintegrable dynamical systems in the separatrix region 
of non-linear resonances. The name results from its ex- 
tensive use in the investigation of chaos, especially the 
mechanisms involved in the transition to global chaos in 
conservative systems. It is derived from the kicked rotor 
model of a one dimensional pendulum and is a Hamil- 
tonian (area preserving) dynamical system. Though the 
map has been extensively analysed Jl^-|l5|l, here we de- 
scribe some relevant details. 

The standard map can be represented by the equations 
of motion: 



K . ,„ s 

Pn+1 =Pn- 7T Sin(27T(7 n ) 

q n +i = q n +Pn+i (mod q = 1), 



(1) 



where K is a real variable which acts as the chaos pa- 
rameter. With unit mass and discrete time, p and q have 
the same dimensions. The phase space of the map is pe- 
riodic in q, by definition, and as a special feature of the 
standard map it is also periodic in p, with the same pe- 
riod 1. K — is the case for a free rotor, a trivial and 
completely regular map. In Fig.l(a)-(c), we show the 
well known sequence of standard map plots for increas- 
ing values of K. The map has several axes of symmetry 
and the unit periodicity in both p and q means that only 
a unit square in phase space need be viewed. The inter- 
val [— i, |) is used for all results in this paper. Fig. 1(a) 
shows the mapping when K = 0.2. Much of the phase 
space is composed of KAM (Kolmogorov-Arnold-Moser) 
tori H stretching horizontally from q = — | to g = | 
and serve to isolate one region of phase space from an- 
other. Periodic orbits in the central resonance can be 
easily seen, as well as a few other dominant non-linear 
resonances (small ellipses) between KAM tori. 

As K is increased, KAM tori that horizontally span 
the phase space are destroyed by resonances and are re- 
placed with smaller KAM islands. Beyond a critical value 
at K » 0.97, the last phase space spanning KAM torus is 
broken and the map becomes globally chaotic (Fig. 1(b)). 
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Now the remaining resonances are clearly visible as sta- 
ble islands in a chaotic sea of trajectories. On reaching 
K = 4, most structure is wiped out (Fig. 1(c)). 



B. KS entropy contours 

A positive Lyapunov exponent, which quantifies the 
exponential divergence in time of two closely neighbour- 
ing phase space trajectories, is a primary definition of 
classical chaos |[(|. For one-dimensional maps such as 
the standard map the positive Lyapunov exponent is 
equal to the Kolmogorov-Sinai (KS) entropy, hxsi which 
is a measure of the rate of information production in the 
system |TJ,[ll|. Thus hxs — only for completely reg- 
ular dynamics. The KS entropy (also called dynamical 
entropy or metric entropy) of a chaotic mapping can be 
calculated using the formula 



J its 



lim 

t — >oo 



t 

E 

n=l 



l0g 2 In, 



(2) 



where l n = y/ (Sp n ) 2 + {Sqn) 2 is the changing distance 
between two initially close neighbouring points, (po,Qo) 
and (po + Spo, <7o + 5qo), in phase space, dp and 5q are 
evolved by iterating a linearised form of the chaotic map. 
This tangent map is rescaled after every iteration as fol- 
lows: the n th iteration of the map produces the values 
Sp n and 8q n from which l n is calculated. These values 
are then rescaled to 5p n = 6p n /l n and 5q n = Sq n /l n 
which are fed back into the tangent map for the next it- 
eration jl2| . Use of the base-2 logarithm in (||) allows the 
entropy to be measured in bits of information. 

The standard map is linearised to give its associate 
tangent map 



Spn- 
5q n ^ 



1 K COs(2irq n ) 

1 1 + Kcos{2iTq n ) 



Sp n 

Sq n 



(3) 



which can be employed in (0) to calculate h^s- 
The tangent map (^) clearly shows that the value of 
has depends on the initial position in phase space {po, qo) 
for the standard map. This is not always the case. 
Iiks is generally used as a global measure of the level 
of chaos in a given system, but this is only valuable if 
all chaotic trajectories in the system can reach into all 
regions of its phase space. Well known examples of such 
systems include the cat and baker's maps [|l9|. How- 
ever for K < 1 the standard map has a predominantly 
mixed phase space in which different chaotic regions are 
not connected. KAM tori act as boundaries so that tra- 
jectories originating in one chaotic region cannot escape 
to another. This isolation inhibits the exponential di- 
vergence of chaotic trajectories so that the positive Lya- 
punov exponent, and consequently hxs, will vary from 
region to region. The kicked top is another example of a 
well known mixed phase space system 20 1. 



To reveal a complete description of the standard map at 
a specific K in terms of KS entropy, many values of hxs 
corresponding to many initial positions in phase space 
can be plotted as a contour map on phase space. This 
has been done in Fig.l(d)-(f). Using (||) and (||), and 
setting t = 10 5 iterations, values for Kks were calculated 
for each point on a 64 x 64 grid spanning the same unit 
of phase space and the same K values as in Fig. 1(a)- 
(c). Shading intensity reflects the relative sizes of the 
KS entropy, hxs — is shown as white on the maps 
while darker and darker shades of grey reflect an increas- 
ing Iiks- The black areas show the largest Kks values 
corresponding to the most chaotic regions of the stan- 
dard map. The resemblance between Fig.l(a)-(c) and 
Figl.(d)-(f) is striking. Stable islands in the classical 
maps translate to stark white patches in the contour 
maps. This is because hxs = for all periods in non- 
chaotic dynamics. The chaotic sea of trajectories in the 
classical maps are also faithfully reproduced as very dark 
patches of similar shape and size in the contour maps. 
All these correlations indicate that hxs presented in this 
way can comprehensively display all the essential features 
of the standard map as it becomes globally chaotic. 



III. QUANTUM KICKED ROTOR 



A. Quantization 

The quantised model of the standard map is governed 
by the Hamiltonian 

H stan (p, q, t) = V - - ^ cos(2tt$) £ S(t - nT). (4) 



This "kicked rotor" describes a free particle of unit 
mass which experiences impulses (kicks) at intervals T. 
Following |H| and p^| , the kinematics are that of finite 
dimensional quantum mechanics with periodic bound- 
ary conditions. Position and momentum space are thus 
discretized, placing the lattice sites at integer values 
<Za = Pa = % for a = — y , . . . , ^ — 1. The dimension 
D of Hilbert space is taken as even and, for consistency 
of units, the quantum scale on phase space is taken to 
be 2irh = j^. Position and momentum basis kets are 
denoted by \q a ) and \p a )- 

Initial states l^o) are assumed to be coherent states 
(minimum- uncertainty states). The fiducial initial co- 
herent state |V , o{oo}) = \qo,Po) is defined as the ground 
state of a special Harper operator 23 1 , which can be dis- 
placed with the appropriate operators to produce all the 
other possible initial coherent states \ipo{ab}) = ka,Pb), 



\i>o{ab}) = ex P exp(-2niap)exp(2nibq)\'tp 0{00} ) 

(5) 
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At time t, the system can be described by the density 
operator p(t) = \ip(t))(ip(t)\ which changes according to 
the evolution equation 

p(t + T) = U sP (t)Ur 1 , (6) 
where the kicked rotor unitary evolution operator, 

Us = exp [-2T ) exp { ***h ) ■ (7) 



B. Von Neumann entropy contours 

Parallelling the classical case, we look for an entropy 
measure to reveal the dynamics of the quantum system. 
Entropy in quantum statistical mechanics is referred to 
as von Neumann (vN) entropy, h V N, (the equivalent mea- 
sure in classical mechanics is the Gibbs entropy) and can 
be defined in terms of the density matrix p of a system 
as 

/i^ = -Tr(plog 2 p). (8) 

h v N is a quantative measure of disorder and can be 
measured in bits. However, the unitarity of Hamiltonian 
dynamical evolution dictates that h V N remain constant 
at all times. The situation can be altered by perturbing 
the system through the interaction with an environment. 
Averaging over the various possible effects of this environ- 
ment will then lead to an entropy increase Ah V N which 
can then be employed to measure the system's chaotic 
nature. This is more than a convienicnt mathematical 
construction. To produce a quantum kicked rotor in an 
experimental situation, the free particle motion must be 
periodically opened up to an environment to allow the 
"kick" to be introduced Q. (Though @ defines the 
evolution operator for free motion experiencing an in- 
stantaneous periodic kick, it is equally valid, as long as 
the free motion is not concurrent 0, for a finite time 
periodic kick which is what is required to realise this ex- 
perimentally.) In doing this the environment itself effects 
the system which naturally causes the entropy increase 
required. 

Thus we choose the environmental coupling to mirror 
the form of the kick in (|4|) viz. the q dependence and in- 
teraction time. We also choose the environment model to 
be a collection of degenerate two-state atoms with a range 
of interaction strengths governed by a normal distribu- 
tion (this is a generalisation of the class of environments 
considered by Schack and Caves so that 

H mt = aCO f^ ] ® J2 & M (n)6(t - nT). (9) 
n=— oo 

Thus during the nth kick the rotor interacts with a 
single two state system with Pauli operator a z (n) and 



interaction strength a. Each of the two-state environ- 
ment systems is equally likely to be in the "up" state 
|t), where <t 2 |T) = IT)j or m the "down" state |!), where 
&z\l) = HI)- Also, a is drawn from a collection of M+l 
independent interaction strengths such that a = ctj for 
j = — 4f , . . . , 0, . . . , 4f . The distribution P aj for aj is 
the normal distribution N(ao, a 2 sd ). 

The combined Hamiltonian for the coupled system and 
environment is thus 

Htot = H st an + Hi n t, (10) 

and the corresponding density operator evolution 
equation is 

P (t + T) = U t ot(a,X) p(t) U t -l(a,X), (11) 

where the combined evolution operator, 

■r-r , ,n ( ~ ia\T cos(2irq)\ » , „s 

Utotia, A) = exp ( ^ qj j U s , (12) 

with A G { — 1,1} is the result of measuring the 
two state environment after each interval to determine 
whether it is an up or down state. (As before, this same 
operator would result if (|p) was turned on for the finite 
time required for an experimental realisation of this sys- 
tem.) The effect of this environmental coupling is to 
produce a multiple stochastic perturbation at the end of 
each time interval. After each interval, there are 2M + 2 
different measurement results leading to 2M + 2 possible 
pure states for the system. Averaging over all these pos- 
sible outcomes in the position basis produces the density 
operator evolution equation 

p xy (t + T)= (x\p(t + T) \y) 

M+l p 

= E -f- £<*i#**(«i>*)7K*)^fo,A)i!/} 

3=1 A=-l,l 

= F(x,y) (x\ UsPWU- 1 \y), (13) 
where 

F{x,y)= E Pa, - (2a jD , sin ^1 sin 

J= i V / 

(14) 

now contains all the perturbation effects due to the en- 
vironment. This causes a vN entropy increase which can 
be determined by tracing over the system so that, 

Ah vN (nT) = -Tr(p(nT) log 2 p{nT)), (15) 

where ~p(nT) is the average density matrix of the sys- 
tem after n time intervals. 

One final step will allow us to see the quantum chaotic 
dynamics. Zurek and Paz [ [To| have conjectured that for 
an open quantum system with minimal dissipation which 
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displays classical chaos, the rate of vN entropy produc- 
tion, h v pf, of its quantum analogue, after an initial de- 
coherence time, td, will rise to a maximum value which 
is solely dependent on the sum of its positive Lyapunov 
exponents. This will continue to be the case until the 
system begins to approach equilibrium when h V N will 
slowly decrease reaching zero at time t eqm . In contrast, 
the entropy production rate of the quantum analogue of 
a regular systems will asymptotically tend to zero well 
before t eqm . Applying this to the standard map, h V N for 
the quantised system interacting with an environment 
should be comparable to the KS entropy of its classical 
(unperturbed) counterpart. Thus for td < nT -C t eqm , 

~ Ah vN (nT)-Ah vN ({n-l)T) 

Kn ~ j, ~ h K s- (16) 

A uniformly spaced 64 x 64 grid of initial coherent 
states (corresponding to an even spread over unit phase 
space) were numerically evolved in time. The maximum 
value of h V N for each evolution was plotted on a contour 
map in a similar fashion to the classical case. Fig.l(g)-(i) 
displays the results for the same three values of K with 
D = 256, a = 0.001, a sd = 0.2a , M = 100 and T = 1. 



IV. DISCUSSION 
A. Quantum-classical correspondence 

There are remarkable similarities between Fig.l(d)-(f) 
and Fig.l(g)-(i). For the same K values, the size and 
location of the various stable islands is analogous, dark 
patches are prevelant in the heavily chaotic regions, the 
axes of symmetry are consistent and the overall complex- 
ity of the dynamics is clearly visible in both. 

There are also differences. The quantum contour maps 
are generally much smoother than their classical counter- 
parts. This is because each initial coherent state in the 
quantum system has a support area causing their evo- 
lution to imitate that of a density of points on phase 
space. Thus neighbouring coherent states will fail to 
achieve dramatically different rates of vN entropy pro- 
duction. Increasing D reduces the supports of the initial 
coherent states as well as reducing the overlap between 
neighbouring states. It was found that this led to a re- 
duction in the smoothness of the quantum contour maps 
making them more greatly resemble the classical contour 
maps. 

The process of calculating these entropy contour 
maps was repeated with variations to N(ao,a.y ar ). For 
aQ,u s d -C K, when any entropy increase is due primar- 
ily to the chaotic dynamics of the system and not the 
interaction, similar results were achieved. 



B. Conclusion 

We have demonstrated that an entropy based approach 
allows classical chaotic dynamics to be accurately mea- 
sured in a trajectory independent way which in turn 
makes it eminently suitable to measure p5|] and anal- 
yse the corresponding quantum chaotic dynamics. We 
have also given numerical support to the Zurek and Paz 
conjecture in a chaotic system which folds phase space, a 
characteristic that their conjecture did not directly take 
into account. 

Entropy measures for diagnosing chaotic dynamics can 
also be employed in other maps. The sawtooth map |26| | 
(which becomes Arnold's cat map jl9| for a specific value 
of the chaos parameter K) does not have a mixed phase 
space so entropy contours would be of little interest. 
However, correspondence can be investigated by compar- 
ing quantum and classical entropy measures for a range 
of K values. Even more interesting is the kicked top EG] 
which is described by a map on the unit sphere. Like the 
standard map, it has a predominantly mixed phase space 
for lower values of its chaos parameter making it ideal for 
comparing classical and quantum entropy contours. The 
kicked top also has a special "order-within-chaos" pjj ] 
feature. In general Hks increases monotonically with 
the chaos parameter K of a given map. However, the 
kicked top has islands of stability reappearing for spe- 
cific higher K values when the map is already globally 
chaotic. This leads to an intricate relationship between 
hies and K which can be compared to numerical results 
for the corresponding Ji v n- We will discuss the results 
for these maps elsewhere. 

Finally, KS entropy is information-theoretically de- 
fined as the rate of production of Shannon entropy (also 
called the Shannon information measure) P8fl . Thus, 
within well defined parameters, we have shown that the 
quantum-classical correspondence of chaotic dynamical 
systems may be realised by viewing the Shannon entropy 
production rate as the classical measure corresponding 
to the quantum measure of the von Neumann entropy 
production rate. These results provide a new diagnostic 
for investigating the chaotic nature of quantum systems. 
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FIG. 1. Quantum-classical correspondence in the standard map. Three values of the chaos parameter K are shown for each 
set of maps clearly showing that correspondence is accurately maintained during the transition to global chaos, (a)-(c) Classical 
map in unit phase space bounded on the interval [—5, \) for both q and p. (d)-(f) Contour plot of KS entropy in unit phase 
space for the classical map. (g)-(i) Contour plot of von Neumann entropy production rate in unit phase space for the quantum 
kicked rotor interacting with an environment. 
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